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Abstract 

We consider the model of self-interacting complex scalar fields with a rigid gauge 
invariance under an arbitrary gauge group G. In order to analyze the phenomenon 
of Bose-Einstein condensation finite temperature and the possibility of a finite 
background charge is included. Different approaches to derive the relevant high- 
temperature behaviour of the theory are presented. 







1 Introduction 



Although Bose-Einstein condensation is well-known to occur for non-relativistic spin-0 
bosons at low temperatures, it is only comparatively recently that the analogous phe- 
nomenon was studied in relativistic quantum field theory |], @ ; S 1> HI- F° r relativistic 
fields Bose-Einstein condensation occurs at high temperatures and can be interpreted in 
terms of spontaneous symmetry breaking. 

The extension and generalization of Bose-Einstein condensation to curved spacetimes, 
and to spacetimes with boundaries has also been the subject of much study. The non- 
relativistic Bose gas in the Einstein static universe was treated by Altaie ||. The gen- 
eralization to relativistic scalar fields was given in Refs. ]7|, |J. An extension to higher 
dimensional spheres was given by Shiraishi ||. More recently a detailed critical exami- 
nation of Bose-Einstein condensation in the static Einstein universe has been performed 



|1C[. Bose-Einstein condensation on hyperbolic manifolds [11], and in the Taub universe 



fl2|l has also been considered. 

By calculating the high temperature expansion of the thermodynamic potential when 
boundaries are present, Ref. ]13[ was able to examine Bose-Einstein condensation in cer- 
tain cases. Later work |I4| showed how to interpret Bose-Einstein condensation in terms 
of symmetry breaking in the manner of the flat spacetime calculations 0,0. A crucial 
feature when the spacetime is curved, or when boundaries are present, is that if the sym- 



metry breaks then the ground state cannot be constant in general [14]]. A later study |[L5| 
showed how interacting scalar fields could be treated. 

We want to continue the considerations on the Bose-Einstein condensation taking as 
a model a set of interacting scalar fields with a rigid gauge invariance under the action 
of an arbitrary gauge group G. In order to examine the phenomenon of Bose-Einstein 
condensation in relativistic quantum field theories, the interesting regime is the high- 
temperature one to which we will concentrate from the beginning. In doing that, different 
approaches may be developed and we will present several of them at different stages of 
the paper. 

For explicity and in order to develop the zeta function approach in combination with 
heat-kernel techniques, we concentrate first on the O(N) self-interacting Bose gas in 
curved spacetime at finite temperature and finite chemical potentials. This continues 
the analysis started in [[[], [|, H from flat to curved spacetime. In section 2 we describe 
in some detail the model and derive briefly the equations necessary for the following 
considerations. The aim of section 3 is to derive the high-temperature behaviour of the 
theory using an expansion around the free field theory. The leading dependence on the 
chemical potentials fi a will be given by the one of the free field. The leading term due to 
interaction may be determined by considering the theory with [i a = (this may be seen 
by dimensional analysis). For that reason, in section 3 we derive the high-temperature 
behaviour of the Bose gas for vanishing chemical potential and for arbitrary N. In that 
case it is no additional complication to allow for nonconstant background field and we 
will consider this general case. Afterwards, in section 4, we concentrate on the 0(2) 
model and we restrict our attention to a constant background field. Because we are only 



1 



interested in the leading dependence of the effective action on the temperature, this is 
actually no restriction because by dimensional grounds it is clear that derivatives of the 
background field will not enter. This has however the advantage that calculations are 
much more explicit in that for example an explicit knowledge of the excitation energies 
is given. We will develop two different approaches. The first one is once more based 
on the zeta function approach, the second one is based on a weak field expansion of 
the excitation energies. In both cases we determine the leading two terms in the high 
temperature expansion. Section 5 contains remarks on the limit N — > oo in the O(N) 
model, especially the influence of higher loop corrections is briefly discussed. In section 6 
we extend our results to general gauge groups. Here we choose the approach based on the 
excitation energies. In that part the chemical potential is treated as a small interaction 
term and we use normal quantum-field perturbation theory. In this approach and for the 
general gauge group we indicate how results containing higher orders in the temperature 
expansion may be obtained and how a nonconstant background field may be treated with 
in a systematic way. 



2 Incorporation of fi in an O(N) model of a self-interacting 
Bose gas 

In this section we will derive the basic equations of a model of N real scalar fields with 
an arbitrary self-interaction coupled to a classical gravitational field. We are especially 
interested in the influence that conserved background charges have on the theory. As has 
been argued in detail [[J, the maximum number of mutually commuting charges is N/2 
for iV even and for definiteness we will always take N even. 

We will consider the O(N) model in an arbitrary ultrastatic spacetime |16| with 



ds 2 = dt 2 - g ij {x)dx t dx j . (2.1) 

x % are local coordinates on the spatial hypersurface £ which is taken to be a .D-dimensional 
Riemannian with a finite volume V. The infinite volume limit can be taken at the end. 
The scalar field action is 

S = Jdt jja x {^(d^ a )(d^ a )-^(m 2 -CR)m 2 -V tree (^)}, (2.2) 

where V tree (Q) is the interaction potential at the tree level. We will assume that this 
potential only depends on $ a through |$| 2 = $ a $ a . Here and in the following the sum- 



mation over a = 1, N is always included. da x = Jg%j{x)d D x is the volume element on 
S. 

Associated with the invariance of ( |2.2j ) under global gauge transformations, there are 
N(N — 1)/ 2 conserved currents J^ b giving rise to the conserved charges Q a b '■ 



Q ab = I da x J° ab (2.3) 
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where 

J° ab = i($ a $ b -<$> b $ a ) . (2.4) 

A quantum state of the system is specified by the eigenvalues of a complete set of mutually 
commuting charges. As mentioned (see @), the maximum number of mutually commuting 
charges is N/2. A convenient set of mutually commuting charges is Q12, Q34, which we 
label by 



Qa Q 



2a-l,2a 



da x f a . 



(2.5) 



The conservation of charges is implemented in the usual way by introducing Lagrange 
multipliers /i Q and defining the generalized Hamiltonian density H by 



H = H — /i Q J° 



(2.6) 



where 7i is the usual Hamiltonian density. Here and in the following, greek indices are 
always summed from 1 to N/2. 

The canonical momenta which follow from (12. 21) are 



and we have 



n 



1 



nji a + -|v$| 2 + -(m 2 + £i?)|$| 2 + W$) • 



(2.7) 



(2- 

Here |V$| 2 



2 2' 

The generalized Hamiltonian density follows simply from 
g^(d^ a )(d^ a ). 

The reason for adopting a Hamiltonian approach is, that this is the easiest way to 
incorporate finite temperature effects in quantum field theory. First of all perform a Wick 
rotation r = it to obtain a (D + l)-dimensional Riemannian spacetime. Since we have 
restricted the spacetime to be ultrastatic, this presents no difficulty The grand partition 
function is expressed as 



Z = J[dU a }[d^ a ] J^dr J^da x {iU a ^ a -n} 



(2.9) 



The path integral in (|2.9| ) extends over all fields $ a periodic in time with period (3 = 1/T. 
Because Ti is quadratic in the momenta, the integration over the momenta in ( |2.9| ) may 
be performed simply by completing the square [ i~7| . This leaves a configuration space 
path integral 



Z= /[d$ a ]exp(-S[$]), 



(2.10) 



where 



J dT da x (<t>2a-l ~ iy a $2a) + ^ (^2a + i^a® 

-^iv$i 2 - \{m 2 + e^)i$i 2 - Vtr,m\ ■ 



(2.11) 
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Rather than deal with the thermodynamic potential, we will consider the finite tempera- 
ture effective action computed using the background field method |18| . Thus we introduce 



some scalar background field $. By expanding about this background in the usual way, 
the one-loop effective action turns out to be 



T = S[$] + ^lndet{A%[<l]} 



s[$] + rw, 



where 



5$ i (a;)5$ i (x / ) 



(2.12) 



(2.13) 



The second term in (|2.12|) arises from the Gaussian functional integration, and rep- 
resents the one-loop quantum correction to the classical action functional. A is a unit of 
length introduced to keep the argument of the logarithm dimensionless. Using (|2.11|) , for 
Sij [4>] we find 



%[$] = (Mij+V^S&x'), 



where 



^ M 2a>2a _ 1 M 2a , 2a J - [ 
the other components are zero, and 



—2ifi a -^: 



2?< j~Jj(% Qq- 

□ - n 2 a + m 2 + £R 



(2.14) 



(2-15) 



VI ($) 



e i e j j 



(2.16) 



Here we introduced e« = $j/|<l ) | 2 . The separation in M and V is a separation in free field 
plus terms arising due to interaction. 

In order to calculate the functional determinant of the operator SV, ^], we will use the 
zeta function definition introduced by Hawking |19[ and by Critchley, Dowker [p0[| . Then 
( [Dp reads 



1 



r = £[$]-- C'(0) + C(0)lnA 2 



(2-17) 



where ((s) is the zeta function associated with the operator SV,[$], equation Q2.14|) . 
We are especially interested in the high temperature behaviour of the theory. In that 



regime heat-kernel techniques have been shown to be a very powerful tool [21, 14] and will 
also be used here. In doing that, one of the encountered problems is the diagonalization 
of [$] [|l"5| , p2]] . For fi a = this may be done |22| and we will first consider this case in 
detail. 
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3 0(iV)-model with vanishing chemical potentials 



In considering vanishing chemical potentials fi a = 0, the problem is enormously simplified, 
because the operator M, ( 2.1 5| ) , is diagonal, My = — + ^R + m 2 ). For that case, the 
eigenvalues of the operator may be represented in the form 



a, j 



3) ■'« 



.)•■ 



with the eigenvalues Xj of the operator 



D = -A + [/($). 



(3.1) 



(3.2) 



For simplicity we introduced here 



UijW = V tJ + ^R + m 2 )S,, 



i.r 



Using a Mellin transformation, the zeta function of Sij[3>] may be written as 

COO = 



r 



bjL?f* < M-[(£) ,, H t } 



— nH\ TrK(x,x,t), 



(3.3) 



where the kernel K(x,x',t) satisfies the equation 



d 

- -^K(x,x ,t) 
limi^(x, x' t) 



(-A + U($))K(x,x',t), 
\g\~' s 8(x, x 1 ). 



(3.4) 



We would like to concentrate on the high-temperature behaviour of the theory. Then 
for n ^ only small values of the parameter t contribute considerably to the integral in 
( |3.3D . So it is reasonable for n ^ to make use of the asymptotic behaviour for t — ► 
of K(x,x',t) ||23|| . As mentioned, in this section we will not assume that the background 
field $ is a constant. However, in order that the expansion that we are going to derive 
is reasonable, <5 has to be slowly varying, together with the condition <C 1, where 

\R\ is the magnitude of a typical curvature of the spacetime. 

The suitable tool to derive the described expansion is the use of the following ansatz 



for the heat-kernel suggested by Parker and Toms ||24 

K(x,x',t) 



a(x, x 



exp 



2t 



-R 



t > x 



A VM (x,x')Q(x,x',t), 



(3.5) 
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where 2a (x, x') is the square of the proper arc length along the geodesic from x' to x and 
Ayw(i, x') is the Van Vleck-Morette determinant. For t — > the function Q(x, x', t) may 
be expanded in an asymptotic series, 



Vl{x, x',t) = ai(x, x')t l 
1=0 



(3.6) 



where the coefficients have to fulfill some recurrence relation. Using the given ansatz, 



75]), it has been shown by Jack and Parker [25], that the dependence of a/, I = 1, oo, 
on the field <§ is only through derivatives of the field. As a result, combining (|3.3|) , (|3~5| ) 
and ( |3.6|) , we will automatically find a derivative expansion in the field $ without any 
effort. 

Although the method in principle works to every order, it gets very cumbersome after 
the leading terms. For that reason we concentrate from the beginning only on the leading 
terms. However, the method will be clear after that. 

First of all one has 



with the zeta function of the spatial section, 

1 r°° 

Ce(x) = -— / dtr^TtKfaxit), 
1 [sj Jo 

resulting from the n = term, and the finite temperature part 

C^) = ^ E r dttS ~ lex v{-(j) n 2 t\TrK(x,x,t), 



(3.7) 



(3.8) 



(3.9) 



where the prime indicates that the n = term of the sum is omitted. The zero temperature 
contribution represented by ( |3.8|) has been analyzed in detail in [22] using an adiabatic 
expansion method, so let us concentrate on the finite temperature part, eq. (|3.9|). Use of 
the heat-kernel expansion ( |3.5|) leads to 



C.00 = E£00 



(3.10) 



i=0 



with 



1 



(4vr)Tr(s) 



Tr < exp 



dtt 



s-l- 



exp 



1 



'2V 

J. 



n 2 t 



.([/($) _ _ R)t 



ad 1 



(3.11) 



We will concentrate only on the first two terms, C°( s ) an d Cz( s ) (Cl( s ) is zero due to 
di = 0), containing the leading terms of the high-temperature expansion. As one realizes, 
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a diagonalization of the potential matrix £/($) is needed. For completeness, the matrix 
S accomplishing the diagonalization is presented in the Appendix A (for details see |22| ). 
Using the results presented there, one first finds 



(200 



E 



dt t 



s-l- 



2vr 



exp 



(47r)-r( S )^ 

exp-f-M 2 *} + (N - 1) exp{-M 2 t} 

where we introduced the effective masses 



j " t>x 



Ml 



1 
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(3.12) 

(3.13) 
(3.14) 



These masses, or more detailed (Mj/3) 2 , play the role of a natural expansion parameter of 
the theory. Performing the integration in ( |3.12| ) leads to Epstein type zeta functions, 



Ef 



E^ + c '- 

i=i 

,r(s + o 



E("l) /i 7^C/ ? (2 S + 20c 2i 



z=o 



lW(s) 



(3.15) 



where the last equality holds for |c 2 | < 1 Cr( s ) is the Riemann zeta function. With 
this, in principle all terms in the high temperature expansion may be calculated. 

However, for simplicity let us present only the case D = 3 and only the three leading 
terms. The relevant results then read 



C°(o) 



(3 



2(4tt) ; 



Jja x [Mt + {N -l)Ml\ 



(3.16) 



and 



c;°(o) 



J^da x 




H " f 







3 fl 



[M 2 + (N - 1)M 



0((M 4 /3) 6 )}. 



For Cz( s ) we continue as indicated in the Appendix A, eq. 
to /r 3 C>((M;/3) 4 ), the result reads 



[Ml + (N — 1)M 2 4 ] 

(3.17) 

. Including terms only up 



v U + k 



(47T 



Vis) 



2tt 



2s+l 



Ck(2s + 1) / da x A 



(3.18) 



with 



N N 



N - 1 
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-□ 



'Ml' 



This gives the contributions 



C 2 (o) 



1 







and 



C 2 (o) 



2(47r)2 v/tt is 



6 



da^A 



7 + In 



47T 



darA. 



(3.19) 



(3.20) 



(3.21) 



(47r)2 V 71 " 

Dealing with a manifold without boundary, the terms in A containing the Laplace- 
Beltrami operator vanish. 

Putting things together, we find for the high-temperature behaviour of the O(N)- 
model the result 

3 



r (i) 



2(4tt 



da x 



45 



N 



2tt 
J 



71 2tt 



[Ml + (N - 1)M 



2^ 



7 + In 



4:71 



+p- 3 0((M i pf)} 



3 /3 

[Mi + (N - 1)M 2 4 + 2A] 
C s (0). 



(3.22) 



For the case N — 1 it reduces to the results found in ^7|, as it should. 

Let us now consider the 7^ 0. Because we will indicate how to develop a 

systematic approach in section 6 for the case of a general gauge group, here only some 
comments. Starting point of the calculation is 



COO 



1 poo 

— - £ / dtf-'Tr {exp[4-[$]t]} 



(3.23) 



with Sij[$] given in ( |2.14j) . As seen in (|2.15|) and ([2.161) , splits into the free case and, 

for small coupling, into small corrections due to the self-interaction. Treating the self- 
interaction using normal quantum-field perturbation theory, the leading term clearly gives 
the theory of N/2 charged scalar fields with associated chemical potentials fi a and the 
effective potential, or up to the order we write it down also the thermodynamic potential, 
for the O(N) model at finite T and finite \x a reads, 



n 



7T 

90 



N/2 



a=l 



(3.24) 
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In section 6 we will generalize this formula from the gauge group O(N) to an arbitrary 
rigid gauge group G. 



4 0(2) model at finite chemical potential 

Let us now concentrate on the 0(2) model with a constant classical background field. For 
N = 2 the analysis is somehow more explicit and systematic approaches may be developed 
at finite background charge. We will present two of them. One is based on an expansion 
around [i = 0, the second one is an expansion in the coupling constant which is assumed 
to be small, this is in powers of the potential V tree and its derivatives, as indicated at 
the end of section III. In both cases dimensional analysis shows, that no other type of 
contributions may appear and that the expansion up to a certain order is complete. In 
this section we will also include the possibility that the manifold has a boundary. 



4.1 Exansion around fi = 



One possibility to treat the 0(2) model in curved spacetime is to follow the approach 
developed for the theory in the Minkowki spacetime 0. In detail the fluctuation operator 
is given by 



S u = 5 4 (x - y) 

S12 = S 4 (x - y) 

S 2 x = 5 4 (x - y) 

S 2 2 = S 4 (x - y) 



-□- fj, 2 +£/"($) -i + 



$ 2 [/'($) $ 2 



$2 



$ $2 



+ [/"($) 

OT 

d 



$=<!> 

$1$ 2 U'{§) $i$ 2 
$2 (|) <J>2 

$1$2 U'($) $l$ 2 



$ $2 



(4.1) 



-□-/i 2 + t/"(<f)-f + 



$ 2 [/'($) $2 



$2 <|> $2 



with [/($) = (1/2) (m 2 + ^)$ 2 + V tree ($). 

Doing the Gaussian integration, the one-loop approximation of the partition sum is 



lnZ« 



_^ tr ln D 2_ 4/i 2 9 2 



V' W + ^-2, 2 )n 



(4.2) 



+ ([/"($) 



Introducing 



7 



(-A + f/"(<5)-/i 2 ) (-A 



(4.3) 
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and 



a 2 = 2 



(4.4) 



this may be rewritten as 

lnZ (1) = -trln U-8 2 T + ad T + <y){-% - ad T + j)\ . (4.5) 
2 

Assuming now a complete basis of eigenf unctions of the Laplace-Beltrami operator —A 

with eigenvalues A n and denoting the resulting quantities of eqs. ([Of), ( |4.4j ), as 7^ and 
a 2 , this results formally in 

InZW = E zero _ point + VI, (4.6) 
where we introduced the infinite zero-point energy, 

p _ E+ ' j + E ~' j (A 7\ 

zero— point / t r> ' \ / 

3 

and the thermodynamic potential 



^"7 



i E [ ln i 1 ~ e~ pE+ -) + In (l - e"^--)] . (4 



.1 



Furthermore, the excitation energies are 



E ±!n = p n + ^f±^f. (4.9) 

The representation eqs. ( f4.7j ), ( |4.8[ ), of the one- loop grand partition sum is completely 
analogous to the case of a free scalar field in curved spacetime, the only change is in the 
excitation energies. Putting U(<&) = (l/2)(m 2 + £-R)$ 2 , this is the free field, the grand 
thermodynamic potential of a free field propagating in curved background is thus found 
directly. 

Also for iV = 2, in general it will be impossible or at least very difficult to treat the 
one-loop contribution, eq. ( fOf ), exactly. However, we are interested in the Bose-Einstein 
condensation and the relevant range for this phenomenon is the high-temperature regime 
to which we will once more concentrate. 

By the mentioned analogy with the free field, one may obtain without any calculation 
the high-temperature expansion in terms of quantities given through the zeta-function of 
the excitation energies, 

C^HE^-r, (4.10) 
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or the associated heat-kernel coefficients, 



E± t 

3 ~ 



\AwtJ 



E bftK 

1=0,1/2,... 



(4.11) 



By introducing E± tZeropoint , fl±, as the quantities in ( |4.7|) and ( f4.8|) resulting from the 
excitation energies i£±j, the results read 



4 ppc±( 4 )+ (4^ 



(4.12) 



; 6|ln/5 + P± + S± 



with 



S. 



± 



E 

r=l/2,l, 



u 2+r 



(3\ 2r (2r)\ 



Air J T(r + 1) 



Cr(1 + 2r) 



(4.13) 



and 




-4+2r 



r(2-r)C«(4-2r). 



(4.14) 



Here, PP(±(— 1/2) denotes the finite part of C±( s ) at s = —1/2. 

Thus we have reduced the task of determining the high-temperature behaviour of the 
theory to an analytical treatment of C±(^)> ec L- Q4.10 ). 

However, due to the very complicated structure of the excitation energies, eq. (|4.9|). 
this is also a very difficult task. Even for the determination of the heat-kernel coefficients, 
eq. ( f4.11| ), it is not clear how to do it for the coefficients b± separately. To explain the 
difficulties, let us write Xj for the eigenvalues of the Laplacian, then we have 



E. 



± 



2fi 2 { 2\j + 



U"(§) - 
U"{$) 



U'{$) 



(4.15) 







l<E>=<3> 





U"(<S>) 



I *=4 



Here the difference to the case of vanishing background charge is seen explicitly. In 
contrast to a second order elliptic differential operator (to which eq. (|4.15 ) reduces for fi = 
0) in this case one has a pseudo-differential operator. Whereas in the first case one could 
use very powerful results of mathematicians and physisists for this kind of operators, to 



our knowledge no results on the heat-kernel expansion, eq. ( 4.11 ), are available. However, 
as we will see in the actually needed sum of Q + and Q- we avoid the mentioned problem 
and the calculation of the coefficient bf reduces to the known case \i = 0. 
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For notational convenience let us write 



U"{$) + 



[/'($) 



U"{$) 



2 



thus the eigenvalues read 



As we see in eq. 



E± d = dj + n 2 ± y/AfiPdj + b. 
we need the asymptotic t — > expansion of 
K(t) = K + {t) + K_{t) 



(4.16) 



/ 1 \ 2 

V4^t 



E W + V)^- 

!=0i,l,... 



(4.17) 



Expansion of the exponential factor containing the square root leads to the cancellation 
of terms containing explicitly a square root, letting us with 



K(t) 



2 EE< 

3 1=0 



t 



21 



(21) 



■(Affdj + b) 1 . 



(4.18) 



The powers in dj may be produced by differentiation of the exponential, which then yields 
a reduction to the heat-kernel connected with dj which is nowadays a very well known 
problem. In doing that, it is seen, that the leading term for t —>■ of the summation index 
I is given by t~2 +l . Thus, reducing our attention to the coefficients bo, b 2 , b\ = bf + bf , 
we find that the following terms may contribute 



with 



A 
A 2 



K(t) = A + A 1 + A 2 + 0(t) 



2 (l - + - [b + /i 4 j t 2 + 0(t 3 )j e ~ 
4/i 2 t 2 (l-/i 2 t + 0(t 2 ))E^e-^, 

3 

4 



(4.19) 



(4.20) 



^¥(i + o(t))E4 e 



2-djt 



Introducing the kernel of dj, 



B(t) = -£e-'' 



(±- 



E ^, 



i=o,§,i,.. 



(4.21) 
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this may be rewritten as 



1 3 

( — J X 2 (c + Ci *3 + [Ci - C /i 2 ]t + [C| - /i 2 Ci ]t 3 



V4vrt/ 

+ [c 2 - /i 2 Cl + -(6 + //)c ]t 2 + Oifi 



A 2 



\4irt 



x (6c /i 2 t + Vcitt + 2[c ly u 2 - 3c /x 4 ]t 2 + 0(^)) , (4.22) 



Using the knowledge of the heat-kernel coefficients Cj, this solves already the task of 
finding the high-temperature behaviour of the considered theory and it can be done in 
principle to any desired order. Summing up the contributions in equation ( 4.22j ) we obtain 
to the considered order 



K(t) 



+2 



{2c + 2cxt^ + 2 [ci + 2c /i 2 



C3 + CijJL* 



£5 +2 



C 2 + -&c 



(4.23) 



Thus the leading terms of the thermodynamic potential read 

r 2 



-5 c » rl -^ c i Cs(3)r3 



12 



Cl + 2coAi 2 T 2 + 0(T) 



For a manifold without boundary this is for the effective potential, 



(4.24) 



V t 



thermal 



45 6 12 



with 



M 



U"($) 



[/'($) 







(4.25) 



For the interaction potential 



we have 



U(*) = i(m 2 + ^)$ 2 + ^$ 4 



,U Z „,-■■■ {i - -) /?+~A$' 



generalizing the flat space result |j, || § to curved spacetimes. 
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4.2 Weak coupling expansion 

Our second approach is based on an expansion in the self-interaction potential. In doing 
that, let us also present another very simple way to determine the excitation energies E±. 
For convenience we introduce here a comlex field $ = 775(^1 +^2)- Then, the action 
( p. 11| ) is given by 



S = J dtj da x {{& +«>$ t )($-i//$) - |V$| 
-(m 2 +eR)|$ 2 |- 



(4.26) 



The field equations for $ and $t follow from S in the usual way by varying with respect 
to $ and $t independently. It is easily seen that 








-<£ + 2z/i$ + /i 2 $ + V 2 $ - (m 2 + £R)<f> - VlJb, 



-$ T + 2i/i$ T + /i 2 $ f + V z $ 



72At 



(4.27) 
(4.28) 



(Here 



tree 9|<J>| 2 

fluctuations about $. We will write 



Vireed^l 2 )-) The excitation energies are obtained by looking at small 
vill write 

* (4.29) 



$ = $ 

and linearize ( |4.27| , p~2"§| ) in ^ . This gives 

= _^ + 2i/i^ + /x 2 * + V 2 ^- (m 2 + £R)^ 



(4.30) 



along with the equation obtained by taking the complex conjugate of (|4.30|) . It is of course 
not possible to solve ( [4.30| ); however, we can expand ^ and in terms of a complete 
set of functions. Let {f N (x.)} be a complete set of eigenfunctions of —V 2 + £i? with 
eigenvalues : 

(-V 2 + £R)f N (x) = a N f N (x) . (4.31) 

We can write 



*(t,x) 
*t(t jX ) 



E^e-^Viv(x), 



A? 



(4.32) 
(4.33) 



A? 



for some independent expansion coefficients A a? and .B/v- Substitution of ( |4.32;Ji~3~3| ) into 
( |4.3C| ) and its complex conjugate leads to the coupled equations 



(E N + /i) 2 - a N - m 2 

-&v;' ree B N , 

(E N - fi) 2 - a N - m 2 - VL 



tree 



vJLM 2 



^tree I ^ I 



.4 



A 



B 



N 



(4.34) 
(4.35) 
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In order for a non-trivial solution for An and we must have = E±n where 



1/2 -, 1/2 



e ±n = {n 2 + m 2 n ± + m\v;> ree ) } ' } 



with 



Ml 



a N + m 2 + V'+m 2 V/' 



tree 



(4.36) 



(4.37) 



This gives the excitation energies of Refs. |J if the result is rewritten in terms of real 
fields and it is the result (|4.9|) 

Once the excitation energies have been determined, the thermodynamic potential fol- 
lows from the usual definition. Alternatively we can deal with the effective potential or 
effective action. The thermodynamic potential Q is given by 



n 



T5>[(1 

N 



-0E+I 



-PE- N 



(4.38) 



(If we deal with the effective potential then there is a contribution from the zero-point 
energy [|], but it can be ignored at high temperature.) If we set V tree = so that the 
theory is free, then it is easy to see from ( |4.36|) and (|4.37|) that 

E± N = (a N + m 2 ) 1/2 ±/j . (4.39) 

The term (crjv + m 2 ) 1 ^ 2 just gives the energy of a single mode /W(x), and is the analogue 
of the single particle energy (k 2 + m 2 ) 1 ^ 2 in flat spacetime. Using ( |4.39| ) in ( 4.38|) results 
in the standard expression for the thermodynamic potential for a free boson gas at finite 
density. 

When Vtree 7^ 0, the evaluation of ( |4.38[ ) becomes more complicated. A general ap- 
proach to this problem can be made using generalized ^-functions, which has been pre- 
sented in the first part of this section. However, it is very simple to obtain the leading 
term in Q coming from Vtree at high temperature. To do this define 



M 2 N 



u 2 N + 5M 2 



where 

K + m 2 ) 1 / 2 , 

vLeim 2 ) + mxum 2 ) ■ 

If we work to first order in the potential it is easy to show that 

1 



5M 2 



E 



±N 



UJ N ± /! + 



2uj 



-5M Z + ■ 



N 



(4.40) 

(4.41) 
(4.42) 

(4.43) 



Using this in (|4.38 ) and expanding to first order in SM 2 results in 



N 

1 



,-/9(wjv+A«) 



(» 



) 



V— 5M 2 j 

_i_ r e ^(^jv+M) _ - 



}+ 



(4.44) 
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The first term on the right hand side of ( f4.44|) is just the thermodynamic potential in the 
absence of interactions, ftf ree , 



free 



N 



1-e 



)■ 



(4.45) 



To evaluate the second term on the right hand side of Q4.44T ), note that 







dm 2 



In ( 1 - e ~^ N±tl) 



JL 

2un 



-1 



This shows that 



d 



- 1 



-1 



+ 



6M 2 



o, 

n o free 



(4.46) 



(4.47) 



The first order correction to the thermodynamic potential due to the interaction can be 
determined from a knowledge of flf ree - 

Expressions for Qf ree can be found in a number of places |^T, jlBfl . The result for 
a manifold with boundary is (for D = 3) 



_rp2 

12 



0i - (m 2 - 2/i 2 



+ 



(4.48) 



at high temperature. 8k are the heat kernel coefficients for trexp 

(-t(-V 2 + £i?)). The 

result for spatial dimensions other than 3 is also easily obtained : 



n 



free 



_2 7r -(- D + 1 )/ 2 r(( J D + 1)/2)Qr{D + 1)T D+1 9 
-7t^ d+1)/2 V(D/2)Cr(D)T d 9 1/2 

-^ d+1 ^ 2 t((D - 1)/2){ R (D - l)T D - 1 



x 



m 2 — 



(D-i)fi 2 ]e } + 



for L> > 3. 

If we consider only D = 3, and use #o = V", then 







dm 2 



and so from ( |4.47| ) 



free 



1 



12 



T 2 V + 



= —T 2 V5M 2 + 



(4.49) 



(4.50) 



(4.51) 
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We therefore have found the high temperature expansion 



n 



--VT 4 
45 



-7T 

2 



-3/2 



+— T 2 (SM : 



12 



m 



Ck(3)T 3 1/2 

2/j 2 )v - e 1 



+ 



(4.52) 



SM 2 is given directly in terms of the interaction potential by (|4.42|) . If we specialize to S 
having no boundary then 0i/ 2 = and #i = — (£ — 1/6)RV. We therefore find 



V, 



thermal 



Q/V 



45 12 



*rC(l^l 



+m 2 + (f - l/6)i2 - 2// 5 



+ 



(4.53) 



in agreement with eq. (|4.25 ). It is possible to calculate higher order terms in the in- 
teraction potential by obtaining the second order correction to the excitation energies. 
However it is clear on dimensional grounds that such terms will only involve lower powers 
of T, and therefore that these terms will be negligible at high temperature. 



5 O(N) model at finite chemical potential in the limit 

N oo 

Let us conclude the considerations on the O(N) model with two comments on the limit 
N -> oo. 

The first comment is on technical simplifications that occur. In the large- iV limit the 
off diagonal elements of the potential matrix V($), eq. ( |2.16| ), may be neglected 0, so 
that the potential reduces to 

Vi^S^^M). (5.1) 



Thus the fluctuation operator splits into N/2 2x2 matrices, 

s Q m = 



S2a-l,2a-l S^a-l^a 



(5.2) 



2 



So it is easily seen, that also the one-loop contribution splits into the sum of N/2 contri- 
butions, 

i N/2 

lnZ (1) = In det (5.3) 
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Thus in the limit of large N, the one-loop contribution is simply the sum of the contri- 
bution of N/2 complex fields with their associated chemical potentials. 

The second comment is on the influence of higher-loop contributions in the limit 
N — > oo. As is well known, at the critical temperature the one-loop approximation is not 
valid any more. This is made obvious by the fact, that the exact one-loop contribution 
would lead to an unacceptable complex critical temperature p8| . The reason is that 
higher loop contribution become significant and all so called daisy diagrams have to be 



summed [28, Let us look in some detail at this summation for the O(N) theory in 
curved spacetime and let us restrict for definiteness on the A$ 4 self-interaction. Then in 
principle, the expansion around the background field leads also to odd powers in <3>, but 
in the limit iV — > oo these are negligible | 28fl . The relevant quantity to tackle with the 



problem is the 1PI self-energy, which one may approximate by [29 

XTN ~ 1 

with the eigenvalues A; of the Laplace-Beltrami operator As. Let us regularize the ex- 
pression by considering 

and lets derive a high-temperature expansion of that expression. Using the method de- 
scribed in section III, one arrives at 

Ca{s) = — — Ca(s) 

XTN 00 / 3 \ 

J2r(s-- + l)a l (A)(2nT) 3 ^ 2l U2s-3 + 2l), (5.6) 



3(4vr)2r(s) 



1=0 



with the zeta function £a and heat-kernel coefficients a/ (A) of the operator As- From 
general zeta function theory it follows, that depending on spacetime properties there might 
exist a pole at s — 1. This divergent piece has to be dealed with in the zero temperature 
renormalization theory. The leading terms of the finite part of the self-energy a defined 
by the finite part PP of C<r( s )> are 

a ^ + ^T ppu(1) (57) 

The leading term is exactly the flat spacetime result, as to be expected for dimensional 
reasons, and so also in curved spacetime the leading orders of the expansion are not 
influenced by summing over all daisy diagrams. 
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6 Extension to general gauge groups 



We have, following Haber and Weldon 0, studied the O(N) model for N even. In this 
section we wish to consider the more general case of an arbitrary gauge group G with a 
set of interacting complex scalar fields with a rigid gauge invariance under the action of 
G. 

Let $ denote the set of complex scalar fields transforming under some <i s -dimensional 
representation of G. We can write 

$ -> U$ (6.1) 

where 

U = exp(i6%) (6.2) 

with Tj the hermitian generators for the Lie algebra of G taken in the <i s -dimensional 
representation for the scalar fields. We will use % — 1, 2, . . . , do where do is the dimension 
of G, and a = 1, 2, . . . , d s . Because we choose to work with complex fields rather than 
real fields here, the action is 

S = J dt J da x {{d^){d^) -m 2 &<$> - £R$t$ _ V tree ] (6.3) 

where $ may be regarded as a column vector of complex scalar fields with d s components. 
Vtree is assumed to be invariant under ( |6.1|) , and for simplicity we will assume that it only 
depends on $ and through the combination |$| 2 = as before. 

A set of conserved currents can be found in the usual way (see Ref. |3(| for example) 
by looking at the change in (^T^) under an infinitesimal gauge transformation. It is easy 
to show that 

Jf = i^Tjd^ - d^&Tj®) (6.4) 
gives the conserved currents, and hence 

Qj = % f daJ&TjQ - $ f T,$) (6.5) 

gives a set of conserved charges. For the case of O(N) these results reduce to those of 
Haber and Weldon [Q] given in Sec. 2. 

A chemical potential /i' can be introduced for each charge Qi, and a Hamiltonian H 
found which incorporates the conserved charges : 

H = H- fi l Qi ■ (6.6) 

It is now possible to integrate out the field momenta in the Hamiltonian path integral as 
before to leave a path integral in configuration space. This results in 

S = J dtj da x {(& + ip?<$>%) ($ - iyiPTfi) - |V$| 2 

+ £R)\$>\ 2 - V tr ee} , (6.7) 



— (m 2 
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as the replacement for (|4.26[) . (If we specialize to 0(N), decompose the fields into real 
and imaginary parts, and perform a Wick rotation to imaginary time, then the earlier 
result in (|2.11 ) is found.) 



In order to calculate the thermodynamic potential we will content ourselves with a 
high temperature expansion. This can be done in two stages. We will first look at the 
case where the chemical potentials all vanish as we did in Sec. 3 for the O(N) model. 
However we will use the technique based on the excitation energies used in Sec. 4.2 for 
the U(l) model. After showing how the thermodynamic potential can be obtained in this 
case, we will then return to the case of non-vanishing chemical potentials and show how 
a perturbative expansion may be used to obtain the leading terms of the thermodynamic 
potential in the high temperature limit. 

6.1 Vanishing chemical potentials 

Setting [i % = in (|6.7|) and varying with respect to $ and $t results in the field equation 

= -$ + V 2 $ - (m 2 + f - y/ ree $ , (6.8) 

and its complex conjugate. As in Sec. 4.2, let $ = $ + ^ and $' = + ^ where $ 
denotes the background field. Linearizing in \1/ results in 

= + V 2 ^- (m 2 + £i?)^- V/ ree (|<l| 2 )^ 

-V^ififli&V + (6 .9) 



and its complex conjugate. By using expansions as in ( f4.32| ) and Q4.33| ), except that A 



1JV 



is a column vector and -Bat is a row vector, we obtain the coupled set of equations 

= [(E%-m 2 -a N -V^ ree )5 aP -Kee^a]A NI3 

-V^9 a 9 p B Np , (6.10) 

= -V£J>\p\A Nfi + B NP [{E 2 N - m 2 - a N - 

-VLe^M ■ (6.11) 

One simple way of obtaining the energy spectrum is by decomposing 

A N(3 = Ajr p + ^fsA ll N (6.12) 

where 

^pAjff, = (6.13) 

and similarly 

B N p = B^ p + %B% (6.14) 

where 

QpBjff, = • (6.15) 
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— aAjj a 


+ $ a {(a- 


|$ 




|$| 





= aB^ a 




|$ 




1*1 



Using these decompositions in ( |6.10| ) and ( |6.11] ) results in 

^", ee Bi] (6-16) 

*| J C4l (6-17) 
where 

a = E 2 N -m 2 -a N -V; ree . (6.18) 

By contracting ( 6.16|) with and ( |6.17| ) with $ Q it is easy to show that a = or 
a = 2| ( l| 2 l / i " ee , where a = occurs with multiplicity d s — 1 and a = 2|$| 2 V/ r ' ee occurs with 
multiplicity 1. The excitation energies are therefore 

E +N = [ m 2 + a N + V; ree + m 2 V;L] 1/2 , (6-19) 
E_ N = [ m 2 + (TN + V; ree \ 1/2 . (6.20) 

These results agree with the /i — > limit of ( |4.36|) and (|4.37|) . 

We therefore find that the thermodynamic potential for \i = is 

= 0) = TJ2 { ^ (1 - e"^) + (d s - l){ln (1 - e^ E - N )} . (6.21) 

It is now possible to make contact with the (^-function method of Sees. 2,3 to obtain the 
high temperature expansion. Instead of pursuing this we will simply obtain the leading 
order terms in the high temperature expansion using the weak-coupling limit described 
in Sec. 4.2. It is straightforward to show that 



Q(H = 0) ~ T^{rf s ln(l-e- /3 ^) + ^-(e^ JV 
N 2^jy 

x[dXee + 2|$|V/; ee ]} (6.22) 

= \{d s +[d s V; ree + m 2 V:U^}^free^ = 0) (6.23) 



where 



n free (fi = 0) = T£ln(l - e'^) 2 (6.24) 

N 

is the thermodynamic potential for a single complex scalar field with no interactions in the 
absence of a chemical potential. In the case of D = 3, with no boundary effects included 

Q free (fi = 0) ~ ~^VT 4 + ^T 2 V[m 2 + (£ - + • ■ ■ . (6.25) 

(It is possible to generalize this to include boundary effects and spacetimes whose dimen- 
sion differs from 4 by using the results quoted in Sec. 4.2.) 
We therefore find 

n(*i = 0) ~ -^d^ + ^T 2 V{d s [m 2 + {i-\)R 

+v; r j + 2|$|V/; ee } + • • ■ (6.26) 
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as the leading terms in the high temperature expansion. For the £7(1) model, which 
corresponds to d s = 2 (although there is only one complex scalar field the dimension 
refers to the real dimension), we recover the \i — > result of (|4.53|) . For the O(N) model 
it agrees with (|3.24j) . 



6.2 Inclusion of the chemical potential 

In order to find the leading term at high temperature when the chemical potential is 
non-zero we will treat the terms in the action (|6.7p which involve the chemical potential 
as an interaction and work to quadratic order in the chemical potential. It is easiest to 
use the imaginary time formalism in which (6/7) is replaced by 

S = ! dt I d<r x {(& + fJ&Tj) ($ - /I'TjQ) + |V$| 2 
J o Jt, 

+ (m 2 + ZR)\<S>\ 2 + V tree } 



(6.27) 



after the Wick rotation to imaginary time. We will define 



S int = £ dt J da x {^(&T^ - - ii j V k ®%T k <&} (6.28) 

as the interaction part of the action. We are therefore perturbing around the fi = result 
for which we know the thermodynamic potential already from Sec. 7.1. 
Up to second order in the interaction we have 



r 



(Sint) 2^ nt ) 



(6.29) 



where terms denoted by (■ ■ •) are to be Wick reduced with only one-particle irreducible 
graphs kept. (See Ref. [P| for example.) The basic result is 



where 



($ a {x)$l{x')) = G af) {x,x') 



■— - V 2 + m 2 + £R \ G a p(x,x') = 8 a p5(x,x') 



(6.30) 



(6.31) 



defines the thermal Green function. (We have dropped the terms in Vf ree here as the first 
order correction due to the self-interaction was considered earlier.) 
We have 



(S int ) = /J, k (TjT k ) af3 I dt I dcr x G a p(x,x) . 



(6.32) 



Because we are only after the leading behaviour at high temperature, we may approximate 
G a p(x,x') by its flat spacetime expression. We have 



G a p(x,x) 



OO 



d u k 



(5 J (27T) 



D 



2im 



+ k +m 
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<5 Q/3 T(47r)- D / 2 (27rT) D - 2 r(l - D/2) 

oo 

+ 2E 



m 



2nT 



D-2 



n=l 



n 2 + 



m V 
27rf) 



D/2-1 ' 



(6.33) 



For high T, taking D — > 3 it is easily seen from ( p.33j ) that 

1 



This leads to 



G a p(x,x) ~ ^T 2 5 al3 . 



(S int )~--TVfi i fiHr(T i T j ) 



(6.34) 



(6.35) 



as the leading contribution at high temperature. 

The Wick reduction of the second term in ( |6.29|) , keeping only the part quadratic in 
the chemical potential, results in 



(6.36) 



( S int) = 4 ^>- 7 (Tj) Q/3 (r j ) 75 dt dt' / da x / da x > 

Jo Jo Jy, Jt, 

d d 

x fflGp-y ( x > xn > dl Gsa ^ x '' x ) 



after a short calculation. If we again use the flat spacetime Green function, with the 
understanding that we will only keep the leading order term, it may be shown from (|6.36| ) 
that 



(Sint) 



-8fx i fj, j tr(T i T j )V(An)- D/2 T(2 - D/2) 

oo 

x E(2vmT) 2 [(27 m T) 2 + m 2 ] D / 2 - 2 . 



(6.37) 



n=l 



Taking the high temperature limit and letting D — > 3 results in 

(SL) * irvwtr^T,) . 

Using (|6.35|) and ( |6.38|) in (|6.29|) and noting that T = /30, we have 

n ~ -^T 2 y/i> i tr(T i T i ) 



(6.38) 



(6.39) 



as the leading order correction to the thermodynamic potential at high temperature when 
the chemical potentials are non-zero. 

It is possible to relate tr(Tjl^) to the quadratic Casimir invariant of the group using 
tr(TjTj) = (dijC 2 (G r) I N)5ij , where dn is the real dimension of the representation Gr of 
the gauge group, and iV is the dimension of the group. In the special case of O(N) the 
result in ( |6.39p reproduces fl3.24| ) given earlier. 
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7 Bose-Einstein condensation 



Let us finally come to the application of the calculations to the phenomenon of Bose- 
Einstein condensation for the case of a constant background field. The high-temperature 
effective action r[$] including the classical part of the background field $ has the form 



£[$] 



+ eR)|<J| 2 - //'<i»T ( //'y; ; <i> + v tree ($) - ^-d s T 



m + U 



6 



90 



(7.1) 



^VV'Tr^T,) 



Thus we find the equations of motion 



T' 2 



m 



$ = o 



(7.2) 



together with its conjugate. Hence the effective action has one minimum with unbroken 
symmetry, $ = 0, and solution(s) with broken symmetry determined by 



det (m 2 +£R- ,rr,,^lj + V; re M + ^ [d s V^) + 2V»Jp) 
For the 0(iV)-model introduced in section 2 this is equivalent to 



0. 



I$l 



A 



Li He 



m 2 + £R + 



N + 2 
72 



AT' 



(7.3) 



(7.4) 



For N = 2 this generalizes the flat space result 



and reduces to the result found in 



The discussion of eq. (|7\4] ) then parallels completely the flat space discussion [|J and 
we will indicate it only briefly for N = 2. At very high temperature one finds the charge 
density 



P 



i dr 

(3 <9/i 



1 



T 2 ji 



(7.5) 



the first piece corresponding to the charge Q in the ground state, the second one, Qi, 
to the excited states. In the symmetric phase, fixing the charge Q and volume V, one 
thus have /i = 3p/T 2 . As temperature decreases, /z will increase and will reach a value at 
which 



3p 
T? 



™ 2 + ^+^ c 2 



Here T c is, corresponding to (|7.4|), the critical temperature 

T2= 18 

c A 



Li 2 {T c )-m 2 -£R 



(7.6) 



(7.7) 
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at which Bose-Einstein condensation and thus symmetry breaking occurs, 
in the ground state one then easily finds 



For the charge 



Qo = Q 



21 



(7- 



which is identical to the form of the flat space result |2|, ^j, [5] . 



8 Conclusions 

In this paper we continued the analysis started in the references [1-5] from flat to curved 
spacetime. Especially, we considered as a model a set of interacting scalar fields with 
a rigid gauge invariance under the action of an arbitrary gauge group G. In order to 
examine the phenomenon of Bose-Einstein condensation we derived the high temperature 
behaviour of the effective action of the theory. To this aim, different approaches were 
developed, one being based on the zeta function approach in combination with heat-kernel 
techniques, the second one being a weak coupling expansion. 

As an application of our results, in the previous section we discussed the Bose-Einstein 
condensation for the case of a constant background field. Qualitatively the properties are 
the same as for flat Minkowski spacetime, however the results ( |7.6|) and (|7.7|) show how 
the critical quantities at the condensation point depend on the curvature. 

The results given in this article mainly include only the leading order correction to the 
thermodynamic potential at high temperature when the chemical potentials are non-zero. 
A more accurate result which includes sub-dominant terms in the temperature can be 



obtained by using the local momentum space expansion of Bunch and Parker [31] for the 



Green function in place of the flat spacetime Green function ( |6.33 ) we have used. It is 



also possible to extend the result to the next order in fx. This will be given elsewhere. 



Acknowledgments 

It is a pleasure to thank E. Elizalde and S.D. Odintsov for discussions and remarks 
that turned into an improvement of the manuscript. K.K. thanks the Department ECM 
of Barcelona University for warm hospitality. K.K. acknowledges financial support from 
the Alexander von Humboldt Foundation (Germany) and by CIRIT (Generalitat de 
Catalunya). D.J.T. would like to thank the Nuffield Foundation for their support. 



25 



A Appendix: Diagonalization of the matrix 



The aim of this appendix is the diagonalization of the matrix £/($). The matrix S which 
diagonalizes the matrix U($) by Udiagi^) = S^ 1 U(^)S is given by 





( 61 


e 2 


e 3 • 


• e N \ 




e 2 


-ei 


• 


• 


s = 


e 3 





-ei • 


• 




\ e N 





■ 


■ -ei j 



(a.i; 



with the inverse 



s- 1 



( 61 


e 2 


e 3 


eiv \ 


e 2 


l-el 
ei 


^e 2 • 

£l 9 




e 3 


^e 3 

ei ° 


l-e§ 
ei 


■ 



1-e 2 
ei / 



(A.2) 



With these results at hand, eqs. ( |3.12| ) and (|3.18| ) for £° and Q are derived. The only 
thing to use furthermore is 



tr < exp 



tr 



-t (-£/($) - -R 

( exp(-tM 1 2 ) 

exp(-tMf) 



tr \ S exp 






-t{U(S)--R 



\ 



v 



o 







S 1 cijS 



(A.3) 



exp(-tMf) J 
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